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ABSTRACT 

Anticipating the construction of the international space station, on-orbit modal identification of space platforms 
through optimally placed accelerometers is an area of recent activity. Unwanted vibrations in the platform could affect the 
results of experiments which are planned. Therefore, it is important that sensors (accelerometers) be strategically placed to 
identify the amount and extent of these unwanted vibrations, and to validate the mathematical models used to predict the loads 
and dynamic response. Due to cost, installation, and data management issues, only a limited number of sensors will be 
available for placement. This work evaluates and compares four representative sensor placement algorithms for modal 
identification. Most of the sensor placement work to date has employed only numerical simulations for comparison. This 
work uses experimental data from a fully-instrumented truss structure which was one of a series of structures designed for 
research in dynamic scale model ground testing of large space structures at NASA Langley Research Center. Results from 
this comparison show that for this cantilevered structure, the algorithm based on Guyan reduction is rated slightly better than 
that based on Effective Independence. 

INTRODUCTION 

Monitoring the vibration characteristics of orbiting space platforms such as the international space station has received 
considerable recent attention. Proposed designs of space platforms that use flexible truss structures to connect operating 
modules and moving elements (i.e„ science and habitation modules, solar panels and radiators) require that activities and 
experiments performed on-orbit not be disturbed by vibrations induced by motors, positioning thrusters, or other activities. 
Therefore, an accurate and complete understanding of the vibration characteristics is desired. With this in mind, the Modal 
Identification Experiment (MIE) was proposed for Space Station Freedom (SSF) to provide frequencies, mode shapes, and 
damping estimates for model verification (Ref. 1,2). 

Before any platform is launched, extensive finite element models will anticipate the structure’s vibration performance. 
In order to verify these finite element models, on-orbit modal identification must be performed. Such testing can not be 
accomplished on the ground, however, due to the platform’s size and flexibility. To address this difficulty, testing of hybrid- 
scale models has been proposed for initial validation efforts (Ref. 3). Even with what verification can be accomplished 
through these means, on-orbit tests will be required for final verification. 

Detecting structural failure due to micrometeroid impact or other failure mechanisms is also a concern. If, for 
example, a portion of the truss structure that will connect the various parts of the Space Station Freedom were to be damaged. 
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extensive astronaut Extra Vehicular Activity (EVA) would be required to find and repair the damaged section. Without the 
ability to locate a damaged section a priori, more risk would be incurred by the EVA astronauts to manually search these large 
structures. One approach which has shown promise for locating and assessing damage in truss structures is to monitor the 
vibration characteristics of the truss (Ref. 4). 

Modal identification of a structure produces knowledge of its vibration characteristics from time history data. This 
means that the mode shapes, the damping, and the natural frequencies are determined. Sensor placement for modal 
identification of large space structures is a relatively recent topic of study with early works presented and published in the 
late 1980’s and early 1990’s. To date most of the evaluation and comparison of proposed techniques to date has been 
accomplished with numerical experiments. One objective of this work is to apply experimental data from a fully-instrumented 
laboratory truss to evaluate and compare a representative set of sensor placement algorithms. The experimental data involves 
nonideal conditions and measurement errors that can not be fully simulated in finite -element-model-based numerical 
experiments. 

This work presents four representative sensor placement methods and four criteria to evaluate each method. 
Frequency response functions (frfs) are used to extract natural frequencies and mode shapes from which displacements were 
used corresponding to the subset of candidate locations selected by each method. The evaluation criteria applied to these 
experimental results are compared with those applied to numerical data from only using the finite element model. Conclusions 
about pretest planning are drawn. Note that this paper presents a portion of a study which provided more detail and also 
evaluated a fifth method (Ref. 5). 

BACKGROUND 

Ideally, a large number of sensors, often accelerometers, would be placed on the platform so mode shape 
characteristics as well as frequencies could then be extracted from time history measurements. However, because the cost 
of accelerometers that will survive the dynamic launch forces and temperature changes in the on-orbit environment is high, 
it is necessary to minimize the number of sensors required to monitor the structure’s modal characteristics. This implies an 
optimization algorithm in which we determine the least number of accelerometers needed and the best places for mounting 
on the structure to accurately describe these modal characteristics. By definition, the normal mode shapes of distinct modes 
of an analytical model are linearly independent. Therefore, one goal in sensor placement is to have the mode shapes extracted 
from the limited number of sensors also be linearly independent. For space platform application, the algorithm must also 
address the characteristically large number of modes with closely-spaced frequencies that these structures possess. 

Sensor number and placement questions arise in other fields. For example, sensor placement applications for HVAC 
systems to control the work environment in large office buildings could address the possibility of the sick building syndrome 
(Ref. 6). Sensor placement problems also occur in monitoring various manufacturing processes at all levels from raw 
materials to finished product (Ref. 7). Sensor placement in buildings is also important for studying the effects of earthquakes. 
This information helps in the design of new buildings and in the reinforcement of existing buildings to withstand the dynamic 
forces that earthquakes produce. 

In each of the applications above, sensors are placed for one of basically three different goals, state estimation, 
optimal control, or modal identification. The third of these objectives, system identification, is the focus of this work, 
specifically the issue of number and placement of sensors to identify specified modes for space platforms. 

All algorithms in this work use a finite element model as a starting point for selecting sensor locations. For pretest 
analysis and test design, this may be the primary available information. Each algorithm attempts a suboptimal approach for 
placing sensors, an optimal solution often only obtainable through an exhaustive search of every different possible combination 
of sensors. However, in most practical cases, an exhaustive search is computationally intractable. The number of possible 
configurations which must be addressed and evaluated in such a search is 

ti- — - — (i) 
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where. 


r\ = number of possible sensor location configurations, 

I = the total number of possible sensor locations, 

T = the number of sensors to be placed. 

In the small problem of 100 possible locations and given ten sensors to place, 1.7 x 10 13 possible configurations result. 

Therefore in this paper, algorithms which optimize defined cost functions determine the best positions for sensor 
placement. In what follows, the sensor configurations are then evaluated using selected criteria to determine the acceptability 
of the results and to compare the merits of various techniques. In ground experiments, sensors are typically mounted 
triaxially, with three perpendicular sensors mounted as one unit. For the purpose of comparison of algorithms, however, 
individual sensor placement is used in this study. 

Table 1 is a summary of sensor placement literature for modal identification. The methods are listed in chronological 
order with the authors in the left-hand column. The methods primarily fall into four categories with exceptions listed in the 
right-hand column. Note that italicized entries indicate where a previously published method was used for comparison with 
the new method presented in the cited paper. In addition to these works which present new techniques, Larson, et. al. (Ref. 
8,9) conducted recent comparative studies for sensor and actuator placement, one using numerical simulations for the same 
cantilevered truss from which we use experimental data in what follows. 

From the methods presented in Table l, one representative method was chosen from each of the first four columns 
for this comparison study. These methods are the Guyan Reduction method described by Penny et al. (Ref. 10), the Effective 
Independence method described by Kammer (Ref. 1 1), the Driving Point Residue method described by Parker et al. (Ref. 12), 
and the Kinetic Energy method described by Chung and Moore (Ref. 13). They are each described in detail in the next 
section. The Effective Independence method has recently become the most widely used method for placing sensors for on- 
orbit applications and the standard for comparison for new methods for identification of large space structures. 

SENSOR PLACEMENT METHODS 

Effective Independence (El) 

Effective Independence is an iterative method which ranks candidate sensor locations according to their contribution 
to the linear independence of the target modal partition. The target modal partition, <|>, is an m x n column matrix of 
eigenvectors for all the modes of interest, where m represents the number of candidate sensor locations matching the degrees 
of freedom (DOFs) of the finite element model and n is the number of "target" modes of interest. The Effective Independence 
method is usually preceded by calculating the kinetic energy contribution of each DOF. The original candidate set of sensor 
locations is then decreased by eliminating sensors which contribute little to the kinetic energy of the target modes. 

Using Kammer’s notation (Ref. 1 1) the Fisher Information Matrix (FIM), A, is computed as <)> T 4>. The eigenvalues, 
X, and eigenvectors, y, are determined for the FIM, 

i4i|/=A.iJf (2) 

Since the columns in $ each represent different modes of the structure, they are linearly independent. Therefore, the 
eigenvalues of A are real and positive. The eigenvectors are orthogonal and represent n directions in an n-dimensional space 
which Kammer calls the Absolute Identification Space. Then, an m x n matrix G is formed in which each column is 
associated with the corresponding eigenvalue of A, 


G=[H]**[4>4r] (3) 

where "**" indicates term by term multiplication (also called the Hadamard or Schur product). An m x n matrix is 
constructed to indicate the fractional contribution of the 1 th sensor location to the j* eigenvalue as 
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Table 1 . Summary of algorithms found for sensor placement for identification. 
Note: methods used for comparison are presented in italics. 


‘Year Author(s) Refer " c ' ' 

Reduction 

Effective 

Driving 

Kinetic 

Other 


Methods 

Independence 

Point Residue 

Energy 

Methods 


78 Shah & Udwadia™ 14 


78 Le Pourhiet & 
Le Letty Ref 15 


'87 Salama, 

et. al Ref 16 


'90 Parker, et. al. Ref 12 


'91 L.M.S™ 17 


'91 Kammer Ref 1 1 


‘92 Flanigan & Botos 8 ' 7 18 


'92 Kammer & Triller R ' 7 19 



‘93 Chung & Moore Ref 


‘93 Yao, et. al.** 7 22 


'94 Liu & Tasker 11 ' 7 23 


‘94 Tasker & Liu Ref 24 


‘94 Brillhart & Kammer Ref 25 


‘94 Hemez & Farhad 26 


'94 Johnson & Mack 1 ** 7 27 


‘94 Pape Ref 28 


‘94 Larson, et. al."* 7 w 


Covariance 
Matrix Norm 


Error 

Sensitivity 




Weighted 

Weighted 

DPR 

KE 



Genetic Alg. 


Perturbation 
Variance & 

DV 



KE& 
Avg. KE 


Chebyshev 

Interpolation 


Eigenvector 

Product 
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where. 


F e = Fractional Eigenvalue Distribution (m x n). 


Finally, a column vector (m x 1) of the effective independence values is calculated as follows, 

‘ t * k 1^ 

F d = 53 F ei j • F E2 j • — : 53 F&) 

1 >=1 ;=1 


( 5 ) 


where, 

E d = Effective Independence Distribution, 

in which the i* term in the column vector is the fractional contribution of the i* sensor to the linear independence of the 
modal partitions. 

Values in the vector E range from 0.0 to 1 .0. A value of 0.0 indicates that the candidate sensor location does not 
provide independent information and is not observable from its location. Conversely, a value of 1.0 indicates that sensor 
location must be retained in order to identify the target modes in the columns of <(). In application, values are between these 
limits, so a threshold must be established to select those sensors contributing the most to the linear independence. The El 
values for all candidate locations are sorted and the location with the lowest value is removed. The process then iterates 
though these steps beginning with determining the new FIM from the finite element model using only the remaining sensor 
candidates. At each iteration only one sensor candidate is removed. 

Guyan Reduction (GU) 

The Guyan Reduction sensor placement method came about as a result of the work of Guyan (Ref. 29), who devised 
a method in 1965 which is used to reduce the size of finite element models, namely stiffness and mass matrices. This 
placement method works by eliminating those DOF for which the inertia forces are negligible compared with the elastic 
forces. This is an iterative process where during each iteration one "slave” DOF is defined as the DOF with the largest ratio 
of stiffness to inertia and the remaining DOFs are termed "master" DOFs. The stiffness and mass matrices for the structure’s 
matrix eigenproblem are reordered and partitioned, 



where, 

[5] = reordered and partitioned eigenvector, 

[S] = stiffness matrix, 

[M] = mass matrix, 

mm = matrix partition corresponding to master DOFs (m x m), 
ss = matrix partition corresponding to slave DOF (s x s), 
ms,sm = vectors m x s and s x m, respectively, 
s = number of slave DOFs = 1, 
co 2 = eigenvalue. 

Henshell and Ong (Ref. 30) present the constraint equation from the slave partition, 

(7) 
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which is used in the strain and kinetic energy equations for the structure to produce the reduced stiffness and mass 


[S*MS m J-[SJlSJ- l [SJ (8) 

where, 

[S*] = the reduced stiffness matrix, and 
[M*] = the reduced mass matrix. 

The new reduced model is ready for the next iteration and the iterations continue until the number of master DOFs remaining 
equals the number of allowed sensor locations. Henshell and Ong (Ref. 30) suggest only the diagonals of the stiffness and 
mass matrices be used to determine which DOFs are slave and master from the assumption that off-diagonal terms in the 
stiffness/inertia ratio matrix can be neglected. 

Penny, Friswell, and Garvey (Ref. 10) investigated three different structures using the Guyan Reduction algorithm 
and concluded that this method selected an appropriate subset when the structure was grounded or had no rigid body modes. 
They also noted that for structures with one or more rigid body modes, the method works adequately but not as well as the 
Driving Point Residue method which they used for comparison purposes. We are generally interested in free/free structures 
with rigid body modes, so we expect this method to perform less adequately on-orbit than our experimental results indicate. 

Driving Point Residue (DPR) 

Another method used to place sensors is by calculating the Driving Point Residue of each DOF. From Parker, Rose, 
and Brown’s paper (Ref. 12), Driving Point Residues are equivalent to modal participation factors, and are a measure of the 
degree of excitation or participation of each mode in the overall response. Driving Point Residues are proportional to the 
magnitudes of the resonance peaks in a driving point frequency response function. The Driving Point Residue is calculated 
for a structure from the following: 


[DPRM< W **[«(«] 


GO) 


where, 

[DPR] = Driving Point Residue values for the structure, 

[cd] = matrix of eigenvalues of target modes. 

Each element of the DPR matrix represents the Driving Point Residue contribution of that DOF in a particular target mode. 
Therefore the i* row represents the contributions of the i* DOF to each of the target modes. Since some of these elements 
can be small or even zero (if at a node of the target mode), the row is averaged and multiplied by the smallest element in 
that row to produce a "weighted” measure of the Driving Point Residue contribution of each DOF. Since a sensor placed at 
a node of a mode would not yield useful information, this biases the method against small or zero elements in the DPR 
matrix. 

Kinetic Energy (KE) 

The Kinetic Energy approach may be considered the most common approach used to place sensors, since a majority 
of "Engineering Judgement" in placing sensors is based on it. Typically, those experienced in experimental procedures will 
look for, among other things, areas of the structure where there is large displacement. These areas would be the best positions 
for the sensors to observe the behavior of the structure with the largest signal to noise ratios. 

Generally, the Kinetic Energy method works by calculating the kinetic energy of each DOF for each mode and then 
using the DOFs with the largest values as positions for sensor placement. The form used herein takes the equation for the 
m x n modal kinetic energy (MKE), 
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[MXEM<M**[M][4>] 


( 11 ) 


and produces a weighted form once again with the weighted average equal to the average value for the row times the 
minimum value in that row. 

EVALUATION CRITERIA 

In order to compare the selected methods, four evaluation criteria are presented. The criteria were designed to judge 
the ability of a sensor set to accurately determine the mode shapes and natural frequencies of the structure. 

Modal Assurance Criteria 

The first two criteria involve the Modal Assurance Criterion (MAC). The MAC is used to present a measure of 
orthogonality between two eigenvectors (or extracted mode shapes), as, 

MAC = ^ (12) 


where, 

<M>j = i* 0 th mode shape vector, 

MAC values range from 0 to 1 with 0 indicating that the two vectors are independent. Generally the MACs are presented 
in a matrix for ij = l,...,n. From this matrix, two criteria are derived. First, the maximum off-diagonal MAC value is 
determined. Second, the root-mean-square (RMS) of the off-diagonal elements is computed. 

Determinant of the Fisher Information Matrix 

The determinant of the Fisher Information Matrix (FIM) was selected as a criteria because it is used as an objective 
function in one method (Ref. 22) and is the inverse of the Covariance Matrix of the estimate error. The estimate error is the 
difference between the real target mode response and the estimated response (Ref. 11,22). The goal is to minimize the 
estimate error by maximizing the determinant of the FIM. The FIM is given by, 

<?-V* (13) 

r 


Where, 

X 2 = variance of stationary Gaussian white noise. 

The same effect is achieved by maximizing the determinant of the numerator of the MAC described above, 

i4=4> r 4> (14) 


SVD Ratio 

A singular value decomposition of the eigenvector matrix, <j), is determined and the ratio of the largest to smallest 
value is computed. The SVD ratio takes into account three quantities: mode orthogonality, the condition of the mode 
expansion problem, and the observability of the modes (Ref. 10). If a singular value is close to zero, the indication is that 
one mode shape is not independent of the others. With increasing independence, the SVD ratio of highest to lowest is 
lowered. 
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NASA 8-BAY TRUSS EXPERIMENTS 


Experimental data from an 8-bay cantilevered laboratory truss was intended to provide a large database for use in 
validation of proposed methods for on-orbit model verification and damage detection in flexible truss structures (Ref. 31). 
The truss was part of the hybrid-scaled structure designed for the Dynamic Scale Model Technology research program at 
NASA Langley Research Center (Ref. 3) and exhibits desired response characteristics of closely-spaced frequencies and light 
damping. Each bay is 0.5 m (19.685 in.) on a side. The truss was mounted to a rigid backstop plate. 

The finite element model for this structure consisted of 32 nodes and 96 DOFs as illustrated in Figure 1. The truss 
was modeled with rod elements that carry only axial load and concentrated masses for the nodes and instrumentation. Node 
numbering starts at the free end. 

Acceleration time history data was collected at each node in each of the three orthogonal reference-axes directions 
using accelerometers and two driving points which were excited with a burst random pattern as detailed in NASA TM-1 07626 
(Ref. 31). In the peaks of the sample frequency response function (frf) shown in Figure 2, representing the response at node 
30 (near the fixed end) in the -Z direction to the excitation at node 4 (at the free end) in the +X direction the first five mode 
frequencies of the structure are seen. The first five mode shapes, extracted using the polyreference method and all 96 sensors, 
and are depicted in Figure 3 with the fixed end of the truss at the top in each case. Table 2 compares the first five natural 
frequencies from the analytical model and experimental results. 


Table 2. Frequency values for the analytical 
model and experimental data. 


Mode 

# 

Analytical 

Freq.(Hz) 

Experimental 

Freq.(Hz) 

Description 

1 

13.925 

13.876 

1” X-Y Bending 

2 

14.441 

14.480 

1 st Y-Z Bending 

3 

46.745 

48.411 

1 st Torsional 

4 

66.007 

64.035 

2 nd X-Y Bending 

5 

71.142 

67.465 

2 nd Y-Z Bending 


Figure 4 shows a representation of the FEA analytical and experimental mode shapes for mode 1. Observe that in 
the analytical model the X and Z DOFs are equal in magnitude, while in the experimental results the X and Z DOFs are 
unequal. This is because the principle bending axes between the two models are not aligned. Also note that the data from 
sensors at DOFs 45 and 84 are suspect. Here 45 can be seen. Small values are predicted by the analytical model in the Y 
direction and these values therefore show the largest percent differences between analytical and experimental results. 

The mode shapes and natural frequencies were first extracted using the full set of frf data. Subsets of these mode 
displacements were used to evaluate each sensor placement method. In response to the question of how the mode extraction 
might vary with the sensor set selected, only the frfs corresponding to a chosen sensor set were also used to extract mode 
shapes and natural frequencies. In general, the percent differences between the modes extracted from the full frf set and those 
from the selected frf set were less than 1 percent for the X and Z DOFs. Because values in the axial direction of the trass 
(Y) are small, larger differences (typically under 6 percent though) were noted between the mode shapes extracted from the 
full set of frfs and the mode shapes extracted from the reduced set. DOF 23 showed a percent difference, however, exceeding 
200 percent, due to a sign difference between the extractions. The natural frequencies extracted using each method were 
identical to two decimal places. 
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Figure 2. Sample frequency response function 
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SENSOR PLACEMENT EVALUATIONS 


Each of the four methods presented previously were used to select a subset of sensor locations from the full 96-DOF 
finite element model. Subsets included 5% (5 sensor locations), 10% (10 sensor locations), and 15% (14 sensor locations) 
of the available DOFs. Ten percent has been suggested as a reasonable number for on-orbit instrumentation (Ref. 32). The 
5% and 15% subsets were presented to show the trends of each method. "One sensor minus 10%" (9 sensor locations) and 
a "one sensor plus 10%" (1 1 sensor locations) subsets were also evaluated but are not presented here due to the generally 
small change in evaluation values. It is important to note, however, that between the 10% and the "one plus 10%" subsets 
for this truss, the SVD ratio and FIM determinant results improved greatly. Figures 5-7 illustrate the sensor locations picked 
by each method for each case grouped by number of sensor locations. For presentation clarity the diagonal members of the 
truss are not shown in these figures. 

Note that those sensors appearing in the 5% subset are in the 10% subset also and similarly in the 15% subset. So, 
as iterations or reductions proceed for these four methods, the next iteration or reduction always yields a subset of the 
previous candidate sensor location set. 

In Figures 5-7 placement patterns can be seen for each method. The patterns of sensor locations chosen by the El 
and the KE methods are similar. They tend to locate sensors at the mid-section and free-end of the truss, which tends to 
maximize the observability of the first five modes. In other words, from the mode shapes presented in Figure 2 the maximum 
displacements are located in these areas. The El method places more importance on the mid-section of the truss than the KE 
method as can be seen by comparing Figures 5 and 6. 

The GU method distributes the sensor locations over the structure more than the first two methods. Note, however 
that it also locates sensors at the mid-section and free-end, just not as closely spaced as those from the first two methods. 
The GU method is also the only method to pick a subset sensor location in the axial (Y) direction from the four methods 
presented. With the differences between the axial displacements from the analytical model and those from the experimental 
results, no adverse effect was noticed in the GU subset evaluation results. 

The DPR method places sensors at the free-end only. This makes the first mode pair indistinguishable from the 
second mode pair. As a result, this method performs poorly. 

To evaluate each selected sensor subset’s performance, the criteria presented above were applied. High values for 
the FIM determinant and low values for the MAC criteria and the SVD ratio produce a better sensor set. An evaluation was 
first performed using all of the DOFs to provide a relative measure of the effect of experimental imperfections for the fully 
instrumented case. These evaluation results are presented in Table 3. Note that the analytical results are slightly better than 
the experimental results due to modeling error and measurement imperfections. Table 3 is presented to demonstrate 
experimental effects in the criteria and to give a standard for results which follow. 


Table 3 Experimental and Analytical evaluation results using all finite element DOFs. 


Evaluation 

Criteria 

Experimental 

Evaluation 

Results 

Analytical 

Evaluation 

Results 

Percent 

Difference 

Max MAC 

6.455 x 10 ’ 

4.982 x 10 3 

30% 

RMS MAC 

2.769 x 10 3 

2.157 x 10 1 

28% 

SVD Ratio 

1.120 

1.073 

4% 

FIM Determinant 

0.985 

0.990 

5% 
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Figure 8. Evaluation results for 5% instrumentation 


Max Mac RMS Mac 



Figure 9. Evaluation results for 10% instrumentation 


Max Mac RMS Mac 



Figure 10. Evaluation results for 15% instrumentation 
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The evaluation results for the 5%, 10%, and 15% instrumentation cases are presented in Figures 8-10, respectively. 
The MAC evaluation criteria show the GU subset to have the most efficient sensor placement followed by the El subset. The 
ICE and the DPR subsets ranked much lower than the first two. It was found in Reference 5 that for grounded structures the 
GU method picks a subset which is generally better than the others and for ffee structures the DPR method tends to perform 
better. Penny et. al. found similar results with the GU method (Ref. 10). Larson (Ref. 9) found that the El method performed 
better than the KE methods studied. 

The SVD ratio changes dramatically between each of the presented subsets. This is partially due to the discrete 
nature of optimization which is employed by each of these methods. Adding or subtracting one or more sensors a large 
change in the identifiability could result. The SVD criteria also selects the GU subset as the best set followed by the El 
method. 


The FIM determinant, however, shows the El method performing better for all subsets. The GU method is shown 
as the next best for the 5% and the 10% instrumentation. The KE subset becomes the second best method in this criteria for 
sensor sets above 10% (beginning at the "10% plus one" subset). A similar result was found in Reference 5 and it was noted 
that generally the best methods picked by the MAC criteria and SVD ratio were also picked as top methods by the FIM 
determinant criteria, but not always in the same order). 

These four criteria were also applied to evaluate sensor placements using a numerical free/free and fixed/free square 
plate (no experimental data available). Results were similar and, as expected, the DPR method was superior to the GU 
method in the free/free case (Ref. 5). 

SUMMARY 

Four representative sensor placement methods for modal identification were evaluated using four evaluation criteria. 
An 8-bay cantilevered laboratory truss was employed to test each method with the evaluation criteria. 

The four methods presented were representative of nineteen found in the literature. From these, the GU method 
produced the most favorable results. The El method performed second best for this structure. It should be noted, that for 
free structures the DPR method does generally better than GU as was found in References 5 and 10. Also the El method 
performed better than the KJE method as was found in Reference 9. The KE and DPR produced evaluation results that were 
similar for 5% and 10% instrumentation. Above 10% instrumentation, the KE method performs better than both the GU and 
DPR methods considering the FIM determinant criteria. 

All methods were designed for modal identification and evaluated for this purpose as well. However, for damage 
detection, different placement approaches and evaluation criteria may be needed. In particular, more axial measurements may 
be required. The axial sensor locations for this truss structure were found to produce relatively large differences between the 
analytical and experimental results. More sensitive instrumentation should be used for these measurements if required for 
damage detection. 
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